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Abstract 

Using multiscale analysis and methods of statistical physics, we show that 
a solution to the iV-atom Liouville Equation can be decomposed via an expan- 
sion in terms of a smallness parameter e, wherein the long scale time behavior 
depends upon a reduced probability density that is a function of slow-evolving 
order parameters. This reduced probability density is shown to satisfy the 
Smoluchowski equation up to 0(e 2 ) for a given range of initial conditions. Fur- 
thermore, under the additional assumption that the nanoparticle momentum 
evolves on a slow time scale, we show that this reduced probability density 
satisfies a Fokker-Planck equation up to 0(e 2 ). This approach has applications 
to a broad range of problems in the nanosciences. 

Keywords: nanosystems, all-atom multiscale analysis (AMA), Gibbs Hypothesis, 
Smoluchowski equations, Fokker-Planck equations 

1 Introduction 

Nanosystems are currently of great interest in the fundamental and applied life sci- 
ences. A major unresolved challenge is to develop a predictive approach to these 
systems that capture the inter-communication among the processes operating on dif- 
fering scales in space and time. The premise of the present work is that one can 
introduce order parameters (slowly varying quantities that capture the essence of 
large-scale bionanosystem phenomena) and then, using Newton's equations for the 
TV-atom system, derive equations for stochastic order parameter dynamics. 



Examples of bionanosystems abound in nature and medicine. Viruses are supra- 
million atom entities with complex structural and functional characteristics, including 
dramatic transitions, interactions with host cells, and self-assembly of subunits. Ri- 
bosomes are of size and complexity similar to viruses, and mediate an important 
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intercellular process - translation of mRNA into proteins. Protein nanoviruses con- 
duct electric currents, allowing some bacteria to exploit oxide mineral grains when 
performing oxidation in the oxygen-poor subsurface. In addition to these natural 
phenomena, scientists are currently developing nanocapsules for the delivery of ther- 
apeutic payloads (such as drugs, siRNA, or genes) to diseased tissues, and for medical 
imaging by equipping nanoparticles with flourescent subunits while other subunits are 
designed to bind with diseased cell membrane-bound proteins. Finally, researchers 
are designing mutated viruses with diminished viralence to serve as vaccines. 

To address all of these applications, one would like to develop predictive mod- 
ules with the ability to efficiently simulate the dynamics of bionanosystems. Such 
programs should include the following characteristics : 

• an underlying all-atom description to evaluate the interaction of bionanosys- 
tems with selected molecules, membranes, or other features in their background 
microenvironment 

• a model that does not require recalibration with each new application 

• an approach that builds in the detailed physical molecular laws and the predic- 
tive power following from them 

• an approach that is computationally feasible. 

Considering this list, we suggest that a multiscale analysis of the equations of 
iV-atom physics will fulfill each of these requirements. Molecular Dynamics (MD) 
is a current state-of-the-art software package that efficiently performs simulations of 
Newton's equations for each of N atoms in a system of interest. An efficient MD 
code, NAMD, has previously been used to simulate a whole virus using a 1024 CPU 
supercomputer, but the process proceeds at a rate of about 1 nanosecond of simulated 
time per day. The typical timescale for a viral structural transition is on the order of a 
millisecond or greater. Thus, the aforementioned MD code and hardware would take 
3, 000 years or longer to attain meaningful results. As bionanosystems evolve due to 
the cross-talk between processes which take place on many scales in both space and 
time, a computational algorithm based on a multiscale approach seems like a natural 
choice. 

The use of multiscale techniques in statistical mechanics beginning with the Liou- 
ville equation has a long history (see [TO], [6], [7], and [H], and more recently [5], [2], 
and [3]). In the present work, we demonstrate several new elements of the analysis. 
First, in our approach, the nanoparticle's internal atomic state, as well as that of the 
microenvironment, are maintained allowing for a more natural, symmetric starting 
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description. Additionally, we utilize a version of the Gibbs postulated equivalence 
of ensemble and long-time averages, following classical results within ergodic theory. 
A precise representation is obtained for the momentum factor in the normalization 
constant for the lowest-order iV-atom probability density in a perturbation expansion 
of a solution to the Liouville equation. As a result, Fokker-Planck and Smoluchowski 
equations are derived which describe the stochastic dynamics of these slow-evolving 
order parameters. These results can then be utilized in the production of an efficient 
software module that can model nanoparticle behavior over long time scales, thereby 
capturing the necessary structural dynamics of a virus. 

In recent investigations (see [2] and [3]), the reduced probability density W was 
shown to obey an unconserved equation of Fokker-Planck type up to 0(e 2 ). The 
derivation of this equation was inconsistent with the mathematical framework of dif- 
ferential equations as the thermal average and the derivatives with respect to order 
d d 

parameters — — and — — do not commute. In the work that follows, we eliminate 

ambiguities regarding the permutation of the thermal average and these derivatives. 
Additionally, the lowest order distribution was previously taken to be independent 
of the conjugate momentum II. This is done in error, causing the lowest order de- 
pendence on the slow variable II to be lost and propagating this throughout the 
multiscale analysis. In Section 3, we rigorously correct these mistakes and establish 
many of the ideas of [3] on a more precise footing by showing that the correction 
to the reduced probability density W indeed satisfies a Fokker-Planck equation in 
conservative form. Prior to this, we show in Section 2 that if the momentum is an 
atomically varying quantity, rather than a slowly varying order parameter, then W 
directly satisfies the Smoluchowski diffusion equation up to 0(e 2 ). In both sections, 
our derivations occur from the starting point of the general kinetic equation so that 
the resulting coarse-grained equations do not arrive from solubility conditions, but 
from a rational expansion of the Liouville equation. 



2 Multiscale Analysis: Smoluchowski Equation 

A central goal of multiscale analysis is to rigorously derive coarse-grained equations 
starting from a more fundamental, final scale theory. The Liouville equation has 
been a common starting point. The challenge is that while the Liouville equation 
preserves probability by construction, it is not guaranteed that a given truncation 
of the equation will be conserving. A re-examination of multiscale analysis for the 
Liouville equation is now carried out to identify potential difficulties of this type that 
may arise, and to set forth techniques to resolve them. In this section we resolve 
probability conservation violations when the momentum is not a slow variable. 

Consider the Liouville equation in a multiscale framework wherein order param- 
eters are introduced. We consider an iV-atom system consisting of a nanoparticle of 
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iV* atoms and a host medium of N — N* atoms. For each atom i = 1,...,N, we write 
Pi, Tj G R 3 , and m ; > as the momentum, position, and mass of atom % respectively. 
In addition, we use the notation T = {r±,pi, tn,Pn}- For each i = 1,...,N, define 
the indicator function 

q _ f 1, if atom i is in the nanoparticle 
[ 0, otherwise. 

For the nanoparticle, we define its total mass 

N 

m = y^m^, 

i=i 

the center of mass 

N 

rrii 



m 



and the total momentum 



N 
i=l 

To begin the multiscale analysis, we first introduce a dimensionless scaling param- 
eter e in the mass terms by writing 

e = — , 3 

m 

where rh is the mass of a typical atom. In the case that all atoms in the nanoparticle 

have the same mass, m; = rh for all i = 1,...,N*, it follows that — = — . Hence, 

m N* 

e f=s (iV*)~ . In this section, we make the following assumptions: 

1. The total nanoparticle momentum does not evolve slowly - P is 0(e°). 

2. The net force on the nanoparticle is not decreased due to cancelation of atomic 
contributions - / is O(e ). 

3. Large migration distances are not a consideration - R is 0(e°). 

As a result, order parameters are 0(e°) and need not be scaled in e, even though 

dR P P 

Newton's equations show that they evolve slowly as — — = — = e— = 0(e). We 

at m m 

note that other scalings would be appropriate to capture different behavioral regimes. 
Let us assume p satisfies the Liouville Equation 



N 



dt 



Pi d „ d 
rrii dri 1 dpi 



p^Cp (4) 
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where we define F{ to be the force on atom i and t to be time. In addition, we 
assume throughout that p decays at infinity (a standard assumption for a probability 
density) so that boundary terms do not appear in the calculations from integration 
by parts. Denote the collection of all atomic positions by T r = {ri, ...,rjv}- Given 
the probability density, p(T,t), we define 



W(R,t) = J A(T*,R)p(T*,t)dT* (5) 

where 

A(T* rl R) = S(R- R*), (6) 
R is the center of mass order parameter, and 

N 

R* = V— rM 



m 

i=l 



is the T*-dependent value of R. Then, using the dependence of W on p, a solution 
to the Liouville equation, we may show that W must satisfy a conserved equation. 
Since p satisfies (BJ, we find 



dW 



^^-^ m ■ fir* ' fin* ) 

/ N 




A ( / p ACR*dr 



Thus, the reduced probability density, W, satisfies 




rrii dr* 1 dp* J OR 



dW d ( f P* 

/ P^dT* . (7) 



dt dR \J m 

Next, we attempt to determine p up to 0(e). The N-atom probability density, 
p(T,t), is then assumed to be expressed as a function of an additional argument, 
T(r, t, •) in such a way that when the last argument is evaluated at R, p is obtained, 
i.e p(T, t) = T(r, t, R). Instead of labeling this new function, we will just extend our 
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previous notation and refer to it as p(T,t,R). This displays the dependence of the 
probability density on multiple scales of motion. Hence, p depends on the all-atom 
descriptive variables T, as well as on R defined by (OQ), the latter an expression of the 
fact that p has indirect dependence on the all-atom state through order parameters 
and thus depends on the all-atom state in several, distinct ways. 

We apply the Liouville operator to p(T, t, R) and invoke the chain rule to find 

®e - -£ - — — ® p 

dt ' rrii dri OR 

i=l 

Using (JTJ) this becomes 



dp r P dp 

— = -C p — . (8) 

ot m oR 

Here, we are writing £ instead of C because these derivatives are taken at constant 

values of R. By introducing ([3D into ([H]), the Liouville equation fll]) transforms into a 

multiscale equation (see [2J and [3] for more details) as 

^ = (£ + eL x ) p (9) 



dt 

N 



where 

4, = -£ 

i=i 

and 



Pi d „ d 

i-1 . |- F i . 

rrii dri dpi 



(10) 



P d 

^ = -^■4- (11) 
m oR 

Again, it must be noted that C and Co, while seemingly exact in definition, differ 
because the differentiation in Co is performed at constant values of order parameters 
R. Additionally, the differentiation in Ci is performed at fixed values of V. Further 
details regarding the all-atom, multiscale analysis (AMA) for the Liouville equation 
can be found in [2], [3], [5], and [H]. The operators ( TT0"]) and (fTTj) differ from that of the 
previously mentioned papers [2J and [3] since in that work the conjugate momentum 
n is treated as an order parameter. In the work that follows in this section, we treat 
this momentum term as a micro-variable instead. 

Assuming the net force on the nanoparticle does not experience cancelation due 
to fluctuating terms, it can be written in terms of the individual atomic forces as 

N 

1=1 

dV 

We let V(r r ) be the iV-atom potential so that — — = — Fi for every i = 1, N. Next, 
we assume that p may be expressed as a power series in e: 

oo 
n=0 
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A set of time variables, defined via t n = e n t } is introduced to capture effects of 
processes occurring on the various timescales. The chain rule implies 

A = yv— . (i4) 
at 2., at n [ > 

We then expand (J9j) using ffT3]) and (IT3|) and separate e scales. Define for n G {0}UN, 

A n = J- - £ n (15) 

where we take C n = for n > 1. The expansion yields the equations 

A po = 0, 

and for n G N, 

n 

A Pn = - ^ KPn-i- 
1=1 

Assuming the statistical state of the system has quasi-equilibrium character, the 
lowest order distribution p is taken to be independent of t . Thus, to lowest order 
we find 

CoPo = 0. (16) 

This implies po is a function of the conserved dynamical variables, notably the total 

d d 

energy H, as well as of R. The latter occurs because the derivatives — — and — — in 

OTi dpi 

Co are to be taken at constant R, and thus CqR = 0. Then, we can define 

iv 



/ 2 \ 



1=1 

and notice that C H = 0. 

Using the entropy maximization principle, one arrives at the nanocanonical solu- 
tion to (1IE]) from [2]: 

Po = = pW, (17) 

where 

Q{P,R) =J A(T* r ,R)e- pH *dT*. (18) 

Here A is defined as in and H* is the Independent value of H. For convenience 
we write t for the collection of slow time variables t = {ti,t 2 , ■■■}■ As is standard in 
multiscale theory, determination of W is delayed until higher orders in the analysis. 
With this, po is seen to factorize into the conditional probability p (i.e. for T given 
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R), multiplied by the reduced probability W for the slowly evolving state of the order 
parameter R. We define the thermal average of a given dynamical variable, A(T) by 



A 



th 



J pA A(T*)dT*. 



(19) 



Now, we will assume that the nanocanonical ensemble obeys the Gibbs hypothesized 
equivalence between the long-time and ensemble averages. More specifically, we utilize 
a classical theorem of Birkhoff ([I] can provide more detail) which states that the 
thermal average of a dynamical variable A(T) and its long-time average are equal. 
Using classical semigroup methods from applied partial differential equations (see 
[8] for more detail), one may show that the linear operator £ is the infinitesimal 
generator of a strongly continuous semigroup on the function space L 2 (T). This 
semigroup is then well-defined and denoted by e c ° to . Hence, in the analysis that 
follows, we will rely extensively on the property : 



lim - 

t— >oo t 



' Cos Ads = A 



tii 



(20) 



for all dynamical variables A{T). Thus, the long-time average or time evolution of a 
variable does not affect the value of its thermal average as defined in ( TT9l) . The survey 
[1] or the classic article [1] can provide more background information and detail from 
an ergodic theory perspective. 

To 0(e) one finds 



A pi = -Aipo- (21) 
Using the previously constructed semigroup e £o *°, equation (f2~T!) admits the solution 

r-t (l 



Pi 



Ai - 



At 



o C 



to 



,C Q (to-t' ) 



Po^o 

\dW P A dW P rr , 



dp 



dt' 



The first order initial condition A\ is, for now, undetermined and has the dependence 
v4i(r, R, £). As a consequence of the cross-level communication inherent to multiscale 
analysis, the behavior of pi at large to provides information about the ^-dependence 
of W, while the analysis of ([7]) provides a necessary condition on A\ that ensures the 
equation determining W is closed. Letting s = t' — to, one obtains 



Pi 



e c oto Ai _ tQ p 







dW 



+ P 



Pf th W - 



dW 
~dR 



to 



r P 

m 



(22) 



In this equation, (3f = — — (In Q), so that / is the force averaged via the nanocanon- 

oR 

ical ensemble. This term occurs because of the dependence of p on R and we will 



S 



dp 

verify the expression for in the Appendix. Thus, using the Gibbs Hypothesis 
(120]) . we find 

f th = lim - / e- CoS fds. (23) 
to^oo t J_ tQ 

Next, we remove secular behavior from the to-dependence in pi, thereby imposing 
the additional condition that p\ remains bounded as to — > oo. Using (122]) . it can be 
seen that if p\ grows in £ 0) it must do so at least linearly. Hence, we may ensure that 

pi does not grow in i by requiring that lim —pi = 0. We then divide by to, take the 

to— >oo to 

limit as to — > oo in equation (122]) . and use fl20|) . Notice that (P) th = 0, as it involves 

r ( p 2 \ 

terms of the form / pi exp — — dp^. Assuming the first order initial data is taken 



in the nullspace of Cq, that is CoA% = 0, use of (l2"0l) yields 

r Pi dW 
lim — = — - — . 

to^oo to Ot\ 

Hence, we find 

dW 

wr° (24) 

and the reduced probability density W is independent of t\. Note that this property 
follows regardless of the choice of Ai in the nullspace of jC . Using this in (122"]). we 
find 



pi = A x + p 



f3f th W 



dW 



o 



e~ Cos ^—ds (25) 



P 



to m 



dR 

concluding the 0(e) analysis, although A\ and W are not yet determined. 



At this point, one would expect to conduct a 0(e 2 ) analysis of the problem and 

dW 

determine an equation for — — . However, this is unnecessary as the correction to the 

ot 2 

reduced probability density depends only on p an d pi up to 0(e 2 ). Instead, define 
for all n = 0,l,2,... 

W n (R,t) = J A(r r ,R)p n (r,t)dT* (26) 
so that, using ((SJ) and (jTBl . we may write 

oo 

W = Y^e n W n . (27) 

n=0 

Hence, we expand W and p in powers of e as in (fIB"]) and (127]) . Using (]T7]) and (TT9|) . 
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the lowest order correction, Wq, can be calculated as 



Uo = / A{T* r ,R) Po {T*,t)dr* 

Ap (T*,t, R*)dT* 

Ap W(R*,t)dT* 
= W(R,t). 

For e — > 0, one may see that W — > W. Hence, as the long time scales tend to zero, 
the correction tends to the reduced probability density. The 0(e) correction can be 
determined using fliZUj) and (ESJ, so that 



W 1 = J A Pl (Y\t)dT* 



A 



dse 



-Cos 



p* 



to 



m 



AAi(T*, R*, t)dT* - / Ap 



dW 

~dR 

o 



f3f th W 



dT* 
o 

-to 



-C s 



p* 



dsdr 



AAi(r*,i2*,t)dT*+ / Ap I e- Cos —dsdT* ) ■ 



p* 



to 



m 
dW 
~dR 



R=R* 



(R,t) + (3f th W(R,t) 



AA 1 (T*,R*,t)dT* + 
AAi(V*, R*, t)dT*. 



to 



-Cos 



P 



m 



th 



ds 



dW 
~dR 



+ (3f th W 



Now, we may write p{T,t) in terms of its expansion up to 0(e). Using (TTT|) and 
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(|25|) in the right side of ([7]), we find 



f P* 

/ A—p(T,t)dT" 
J m 



[ A—p(T*,t,R*)dr 
J m 

f P* 

/ A— (po + epi) dT 
J m 



Ap^W(R*,t)dT* 
m 



A— Ai(r*, R*, t)dT* 
m 



+e 



JD* 

Ap^ 

m 



p* 



t 



dW 



- Cos —ds Pf th W{R*,t) - —{R*,t) dT* 



dR 



(^) W(R,t) + e [ A^A 1 (T*, R*, t)c?r* 
\m J J m 



p* 



m 



p* 



+e — ■ / e- CoS —ds 



to 



m 



tii 



e I rA^i(r,F,t) + e- 7 



rn 



(3f h W(R,t)-^(R,t) 
(3f th W(R,t)-^(R,t) 



where the diffusion coefficient is 



7 



(P(0) • P(s)) th ds 



(28) 



-to 



and we use the notation P(s) = e C ° S P. Thus, ([7]) becomes 



dW 
~dt 



—e 



d_ 
dR 



A— Ai(r*, R*,t)dT* ) fe 
m 



,d_ 
dR 



7 fdW 
~dR 



- f3f th W 



Using the^expressions for W and Wi,we can expand the reduced probability density 
as W = Wo + eW\. Then, isolating Wk terms and imposing the condition that A\ 
must stay bounded for large R, we find that this equation is closed only if A\ — 0. 
Thus, up to 0(e 2 ), the conserved equation (JTj) becomes the Smoluchowski equation: 



dW 
~~dt 



dR 



7 







rh? V dR ^ 



(29) 



Hence, in the case of fast-evolving momentum, due to P and / being 0(e°), the re- 
sulting behavior of the reduced probability density is governed by the Smoluchowski 
equation (|29|) . In the next section, we alter these assumptions on the behavior of 
nanoparticle momentum and determine the corresponding changes in the structure 
of the equation for W. 



3 Multiscale Analysis: Fokker-Planck Equation 

In this section, we again use multiscale techniques to show that under similar circum- 
stances, the correction to the reduced probability density, W, satisfies a Fokker-Planck 
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equation. In this situation, the momentum is not considered an atomistic variable, 
but instead as an order parameter. Hence, we define T, m, R, and P as before, but 
reformulate the problem to allow for the differing behavior of this slowly-evolving 
quantity. 

To begin the multiscale analysis, we again introduce a dimensionless scaling pa- 
rameter e in the mass terms by writing 

m , , 

e = - 30 
m 

where m is the mass of a typical atom. In this section, however, the assumptions 
on the system of interest change. We are now interested in significant migration dis- 
tances on the order of the nanoparticle diameter, which we take to be 0{e~^) ) and 
hence scale R to be 0(e~z). Additionally, under the assumption that the system 

is near equilibrium, the nanoparticle kinetic energy, — is 0(e°). Using the mass 

2m 

ratio scaling, this implies that P = 0(e~i), as well. Finally, we assume that the net 
force on the nanoparticle is reduced due to cancelation of atomic contributions, thus 
causing the momentum to evolve slowly. Hence, / is assumed to be O(ea). A more 
detailed description of these assumptions can be found in [3]. 

Under these considerations, define the scaled order parameters $ and n by 

$ = e^ (31) 

and 

n = e^P (32) 

respectively. The scaled net force / can then be written in terms of the individual 
atomic forces as 

N 

/ = e-iJ]FA (33) 
dV 

and we let V(T r ) be the iV-atom potential so that — — = — Fj for every i = 1, N. 

Let us assume p satisfies the Liouville Equation (j4]) where we again consider Fi 
to be the force on atom i and t to be time. In addition, we denote the collection of 
all atomic positions by T r = {r±, r^}. Given the probability density, p(T,t), we 
define 

H/($,n,t) = J A(r*^,U)p(V\t)dT*. (34) 

A(r*,$,n) = 5($-$*)5(n-n*), (35) 

N 

1 \ ^ rrij „ „ 



where 

and the terms 



m 
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and 

JV 



i=l 



are the Independent values of $ and EL Then, using the dependence of W on p, 
a solution to the Liouville equation, we may show that W must satisfy a conserved 
equation similar to that of the previous section. Since p satisfies (jl]), we find 



dW 
~~dt 



J A^(r,t)dr 



N 



,i=l 

JV 



mj dr* dp* 



P 1 m . Q r * ^ gp* ) ^ 



J=i 

JV 



P* OR* ^ dR*\ i dA ^ 
^— ' nii or* dp* I <9$ 



,i=i 



p* OP* _ dP*\ i dA 

£ • ITT + F i • ITT) e ~ 2 7^T dT 
/ —f nii dr* dp* J all 

£3 A ( / p ACR*dT* ) +e^A f / pA£P*dT< 



d<S> \J r J dU 

i d { f . P* A i a 7 '• 



Thus, the reduced probability density, VF, satisfies 



i=l 



"" ^r/^rfrA-^l /M/rfT-l (36) 



at a$ \j m J du 

which is in conservative form 



Next, we conduct a multiscale analysis in order to determine p up to 0(e). Sim- 
ilar to the previous section, the N-atom probability density, p(T,t), is assumed to 
be expressed as a function of two additional arguments, T(T, t, -,-) in such a way 
that when the last two arguments are evaluated at $ and II, p is obtained, i.e 
p(T,t) = T(r,i, $,II). Instead of labeling this new function, we will just extend 
our previous notation and refer to it as p(T, t, $, II). 

We apply the Liouville operator to p(T, t, II) and invoke the chain rule to find 

dp _ _ pi_ ^_^P__y^ F , ^ n dp 

dt ' uii dr; <9$ ^ % dp, dU 

i=l i=l 
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Using (13~TT) and (|32|) this becomes 

°P n lP dp r dp .„_. 
oi m o<±> cm 

Here, we are writing Co instead of C because these derivatives are taken at constant 
values of $ and EL By introducing f )30|) and fl32|) into fl37|) . the Liouville equation (jl]) 
transforms into a multiscale equation as 

^ = (C + eC 1 )p (38) 



where 

£ = - y \— ■ — + Fi ■ — 

and 



(39) 



Again, it must be noted that C and Cq, while seemingly exact in definition, differ 
because the differentiation in £ is performed at constant values of order parameters 
$ and EL Additionally, the differentiation in L\ is performed at fixed values of E. 
Unlike the previous section, the operators ( flOl) and ( flTI) are now the same as that 
of the previously mentioned papers [2] and [3] since the conjugate momentum II is 
formulated as an order parameter. 

Next, we assume that p may be expressed as a power series in e: 

oo 

P = J2pne n (41) 

n=0 

A set of time variables, defined via t n = e n t, is introduced to capture effects of 
processes occurring on the various timescales. The chain rule implies 

!=£«•!-• (42) 

We then expand (158"]) using (T4T]) and ( 1421) and separate e scales. Define for n e {0}UN, 

A n = J- - £ n (43) 
ot n 

where we take C n — for n > 1. The expansion yields the equations 

A Po = 0, 

and for nGN, 

n 
i=l 
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Assuming the statistical state of the system has quasi-equilibrium character, the 
lowest order distribution p is taken to be independent of t . Thus, to lowest order 
we find 

CoPo = 0. (44) 

This implies po is a function of the conserved dynamical variables, notably the total 

d 

energy H, as well as of $ and II. The latter occurs because the derivatives — — and 
d 

— — in Cq are to be taken at constant $ and II, and thus £o^ = A)!! = 0. As before, 

OPi 

we define 

JV 



i=i v ' 



and notice that C H = 0. 



Using the entropy maximization principle and proceeding as in [2], one arrives at 
the nanocanonical solution to 



Po = ^ = pW, (45) 

where 

Q(/3,$,n) = J A(r*,$,n)e-P H *dr*. (46) 

This form of the nanocanonical solution is slightly different from that of [2] since it 
was stated in that article that Q is independent of II. We find that this is not the 
case and determine the exact manner in which the IT dependence can be computed 
in the Appendix. Here, A is defined as in fl35|) and H* is the Independent value of 
H. We label t as the collection of slow time variables t = {£1,^2, •••}• As before, we 
define the thermal average of a given dynamical variable, A(T) by 

A th = j * Ap A(T*)dT* (47) 

and use the Gibbs Hypothesis: 

1 f° 

lim - / e- CoS Ads = A th (48) 

t— >oo t J _ t 

for all dynamical variables A{T). Notice that the thermal average, and thus the state- 
ments fj4T|) and fj48l . depend upon the new order parameter n because of fl35|) . 

To 0(e) one finds 

A pi = -A1P0 (49) 
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and using the semigroup e c ° to defined in Section 2, this equation admits the solution 

: e Coto li - 



Pi 



A 1 



to 



e Co(to-t' ) Aif)odt r o 

dW 



e c ° t °A 1 



to 



,Co(to-tf ) 



,Co(t -t' ) 



n ap 

aw n 



p- 



n jsw 

rh P 5$ 

n 



J an j p an 



dt'n 



where we have used the results 



and 



- -(3pf n W + -p— 
oti th m a<P 



avy n an/ 

+ /3/p- W + //)— - 

m an 



dp 

an 



/3p 



n 



These derivatives will be verified in the Appendix. The first order initial condition A\ 
is, for now, undetermined and has the dependence Ai(T, H,t). As a consequence 
of the cross-level communication inherent to multiscale analysis, the behavior of pi 
at large to provides information about the ^-dependence of W, while the analysis jof 
provides a necessary condition on A\ that ensures the equation determining W 



is closed. Letting s = t' — t Q , one obtains 



Pl = e^Ar-top —-tf-W- / e~ Cos (f - f th ) ds 

oti m J_ tn 



~t 

n dW dW f° 

-^-W_ to e /ds - (50) 



Next, we remove secular behavior from the to-dependence in p%. As before, we 
assume that pi remains bounded as to _ > oo. Using (1501 . it can be seen that if pi 
grows in to, it must do so at least linearly. Hence, we may ensure that p\ does not 

grow in to by requiring that lim — p\ = 0. We then divide by to ; take the limit as 

to-s-co to 

to — y oo in the above equation, and use (HHl) . Assuming the first order initial data is 
in the nullspace of C , this yields 

Pi dW n dW , h dW 
hm — = — • — — - / 



Hence, we find 



io^-oo to oti rh a$ an 



Al - m 1 + m a$ + ; an " u (5ij 



and W satisfies a Liouville equation in (t l5 $,n) space. Note that (IBTj) follows re- 
gardless of the choice of Ai in the nullspace of Co- Using this in (1501 . we find 
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Px = A x -p 



n dw 

(3-W + — 
m all 



f e- c ° s (f - f th ) ds 

J -to 



(52) 



concluding the 0(e) analysis, although A x and W are not yet determined. 



As before, the correction to the reduced probability density depends only on p 
and pi up to 0(e 2 ). Hence, define for every n — 0,1, 2, ... 

(53) 



(54) 



w n ($,n,t) = J A(r*,$,n) Pn (T*,t)dr* 

so that, using (134")) and (j4Tj) . we may write 

oo 



n=0 



In addition, we may expand W and p in powers of e as in (HT1) and Using 
the lowest order correction, Wo, can be calculated as 



Wo = J a (r* , $ , n) p (r* , t) dr* 

= J Ap (r*,t,$*,n*)rfr* 

= J Ap w($*,n*,t)dr* 
= w($,n,t). 

For e — > 0, one may see that W — >• W. Hence, as the long time scales tend to zero, 
the correction tends to the reduced probability density. The 0(e) correction can be 
determined using ( 147|) . Notice that (jUJ) implies {e~ Ca ' r A) = A th for any finite value 
of r 6 i. Using this and (|52|) . we find 



Wt 



A Pl (r*,t)dT* 
A 



o 



AAi(r*,$*,n*,t)dr* 

- Api — — + P—W 



8U 771 

AZi(r*, $*, n*, t)rfr* 
AAi(r*, $*, n*, t)rfr* 
AAi(r*,$*,n*,t)dr* 



) ■ |° e- £oS (/ - /») dsdV 



(*,n)=(**,n*) 



(f(s)-f th ) th ds 



to 



dw n 



+ P—w 



dU fh 
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where we use the notation f(s) = e~ CoS f. Now, we may write p(T, t) in terms of its 
expansion up to 0(e). Using (H5j) and (1521) in the first term on the right side of (136]) . 
we find 



A—p(r*,t)dT* 
m 



A—p(r*,t,<$>*,u*)dr* 

m 

A— (p + epi) ^* 
m 

Ap— m$*.n*,t)dr* + e / a— Ai(r*,$*,n*,t)dr* 

m } m 



—e 



Ap^ /5— W + — - 
m \ m oil 



5v^($,n,t) + e5 / AAi(r*,$*,n*,t)dr* 

m m / 



n / o rr cW 

m \ m oil 



(f( S ) - nds 



to 



m m J 

Similarly, the second term on the right side of ([7]) becomes 

Afp(T*,t)dT* = I Afp(T*,t,^*,U*)dT* 



J Af(p + e Pl )dT* 

J Ap/w / ($*,n*,t)dr* + e J A/A!(r*,$*,n*,t)dr* 

dse~ c ° s (f - f th ) dT 



IT* dW 
Apf[(3-W + ^- 
m oil 



to 



f th W($,Tl,t) + e / AfA^T*, $*,!!*, t)dT* 



-e [{3-W+— - 
m oil 



o \ th 

th\ 



/(O) • {f (a) - f th )ds 



to 



— n dw 

f m W(^,U,t)+e I AfA 1 (r*,<5>*,Il*,t)dr*-e6[(3--W+ — 



where 

Thus, ( 1361) becomes 
dW d /n 



[(/(O) • f(s)) th - / tfe " f h ] da. 



(55) 



9* 



_ e J_. rim n, i))+ ^A 



- / A/Ai(r*, $*, n*, twr* + 6> w + 

m on 



aw 
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Using the_expressions for Wq and W\,we can expand the reduced probability density 
as W = W + eW\. Then, isolating Wk terms and imposing the condition that A\ 
must stay bounded for large <& and IT, we find that this equation is closed only if 
A\ = 0. Thus, up to 0(e 2 ), the conserved equation (1361) becomes the Fokker-Planck 
equation: 



dW 



d_ 



n 



-w 



771 



_d_ 

m 



f th W) +e 



,d_ 

an 



n d 

m all 



W 



(56) 



Hence, under the assumption that momenta have (comparatively) large values and 
evolve slowly, i.e., P = 0(e~i) and / = O(ea), the reduced probability density obeys 
a Fokker-Planck equation given by ( |56l) . 



4 Appendix 

We first verify the $ derivative of p used in the derivation of both equations. In 
Section 2, the variable R is used instead of $ (as order parameters are unsealed with 
respect to e), but the statements that follow can be applied exactly to R in the same 
manner as $. We claim 

-0pf h . 



<9$ 

Using (ITTj) or (H5|) . we see that 

dp _e~ fiH dQ 

d 

Notice that — ( e~ pH ] = since $ derivatives are to be taken at constant values of 
T. Using f fl8|) or fj46|) . we see that 



a$ ./ a* 

OA 



-e 



J dR* v ; 



= e 2 

Now, since this integration is performed over all values of P*, and hence T* is not fixed, 

the energy depends upon the center of mass through the potential function. Notice, we 
dV 

may calculate — — as V depends upon R implicitly in the following manner. For every 
3R 

i = 1, N, write the residual displacement of each atomic position as = <Xj + RQi. 
Here the variables {R, ax, &n-i} constitute a complete set of variables as we 
may write in terms of each of the other Cj using the constraint 

N 



^(7^ = 0, 



i=l 
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which follows by the definition of the residual coordinates. This change of variables is 
just a set of linear functions with constant coefficients, hence the Jacobian is constant 
and 

dV ^dV a i 

i=l 

Since the kinetic energy is independent of T r , we find 

dR* K 1 H dR* F J 
Using this, the derivative of Q becomes 



dQ 
<9$ 



e~3 J Ae~P H * Pe^fdT* 



= J Ae- 0H *f*dT*. 

Hence, using the modified Gibbs Hypotheses in either section ( 120]) or (}15]) . we find 

dp e-? H dQ 



<9$ Q 2 <9$ 

1 dQ 



Q d$ 
-pP J ^Ae-^f*dr 

-pp ( [ Apfdr* 



= -Ppf th 

and the verification of this derivative is complete. 

Next, we verify the n derivative of p used in Section 3. We claim 

dp n 

all m 

First, we use (j46p and rewrite Q in terms of its n dependence using the Fourier 
transform. Since the potential and kinetic energies depend exclusively on ri and 
Pi variables respectively, we separate Q into two different integrals involving these 
variables. In addition, we separate the momentum integrals into the momenta of 
particles in the host medium, denoted T = {pi : Q% = 0} and those in the nanoparticle, 
denoted Fx — {pi : 0i — 1}. Using the notation 



N 



k = y*-{i 



2m,- 

i=i % 
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and 

JV 2 
i=l 

for the kinetic energies of the host medium and nanoparticle, respectively, we write 
Q as 



W,$,n) = 




Keeping the T* and Tq integrals as they are (notice further that the Tq integral is 
constant), we focus on the T\ integral. We may again split the Tl integral into one 
each in the x, y, and z directions. We will consider the integral in the x-direction, 
labeled I x , and state that the results we obtain will follow for the integrals in the 
other directions in the same manner. 



Now, we relabel the momenta in the nanoparticle, 1 through M, and write their 
re-coordinates as X\ through %. Then, write the x-directional Dirac mass 5(H X —II*) 
using the Fourier transform of the function f(k) = 1 as 



5(u x - n* x ) 



cxp 



M 



ik j IL X - y/e ^2 x i 

i=i 



dk. 



V27T J-oo 

Then, using the Inverse Fourier Transform on the resulting Gaussians 

1 r°° 1 ' ~ 2 



exp(— Xk 2 ) 

and (j5D|) . we find 
1 



, exp (—ikz) exp I ] dz. 



2tt 
exp 




cxp 



M 



ik | U x - y/e ^ xi 

i=i 



dk 



P > 

-x-, I • • • exp 



2m 1 1 



exp (ikH a 



2m M 



" 4, 1 



'27T -oo 

y exp (— iky/ex exp 



/exp (— zfc\/exi) exp ( — x\ ) dxi 
V 2mi / 



^ x 2 I,/, 



2r«M 



JkTl x 




2(3 



m M 





2m 1 
dk 
evriMk 2 ' 



dk 



2^ 
C\ exp 



m M 



2rh 



2tt{3 y_ 

n 



exp (ikH x ) exp 



efc 2 
2^ 
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where 

/ mi j m M / 27T/3 

1 ~ V 2tt/3 " ' V 2tt/3 ' V m ' 

We extend this in the y and z directions and multiply to find 



J 8(U - U*)e-^ K *dTl = (Ci) 3 exp (~^n : 



Thus, we can write 



where 



and 



Q(/3,$,IT) = Q 1 ($)(C 1 ) 3 C 2 exp ( -^-II 5 



c 2 = / e-^dr ■ . 



Finally, p can be expressed in the form 



Taking a II derivative in (|58|) . which must be done at fixed values of T, we find 

dp n 

oil m 

and the verification of this derivative is complete. 
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